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Abstract 

A thermodynamic analysis of the black hole solutions coming from 
the Einstein-Maxwell-Dilaton theory (EMD) in 4D is done. By consider 
the canonical and grand-canonical ensemble, we apply standard method 
as well as a recent method known as Geometrothermodynamics (GTD). 
We are particularly interested in the characteristics of the so called phan- 
tom black hole solutions. We will analyze the thermodynamics of these 
solutions, the points of phase transition and their extremal limit. Also 
the thermodynamic stability is analyzed. We obtain a mismatch of the 
between the results of the GTD method when compared with the ones 
obtained by the specific heat, revealing a weakness of the method, as well 
as possible limitations of its applicability to very pathological thermody- 
namic systems. We also found that normal and phantom solutions are 
locally and globally unstable, unless for certain values of the coupled con- 
stant of the EMD action. We also shown that the anti-Reissner-Nordstrom 
solution does not posses extremal limit nor phase transition points, con- 
trary to the Reissner-Nordstrom case. 

Pacs numbers: 04.70.-s; 04.20.Jb; 04.70.Dy. 

1 Introduction 

Since the discovery made by Hawking [ ] of the thermodynamic properties 
of the black holes related to the quantum phenomena, many were the interests 
on studying the properties of various kinds of solutions obtained from General 
Relativity and their modifications. In close analogy to the usual thermodynam- 
ics, the black hole thermodynamics is based upon four basic principles, the zero 
law and the other three laws, which are all analogous to the usual ones [2, 3]. 
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Therefore, a new black hole solution can be physically interpreted also by the 
analysis of their thermodynamic properties. Furthermore, we can analyse the 
thermodynamic stability of a new solution through its properties. 

There are several methods to study the thermodynamics properties and sta- 
bility of a black hole. A basic reference on this is the work of Davies [ ]. We 
can mention other more recent methods, such as the Geometrothermodynamics 
[l ] , and the Hamiltonian thermodynamics [ ] . For the present work the hrst 
two methods will developed, always in parallel so we can compare their results. 
We want to stress that the main goal of this paper is to analyze the thermody- 
namics properties and stability of the solutions known as phantom black holes, 
specifically those coming from the Einstein-(anti)Maxwell-(anti)Dilaton (EMD) 
theory. These solutions come from the minimal coupling of the Einstein-Hilbert 
action with a scalar field, that could be either dilatonic or phantom, which at 
the same time this is coupled with Maxwell field, that can be a spin-1 normal or 
phantom field. The phantom term furnishes the contribution of negative energy 
density, which justify the nomenclature. 

In order to develop the analysis of the thermodynamic properties of this class 
of black holes, let us briefly illustrate the interest to studying phantom solutions 
in black hole physics. The programs of evolution of our universe, specially the 
ones for the spectrum of anisotropics of the cosmic background radiation on the 
one hand and for the relation magnitude versus red-shift of the supernovae type 
la on the other, have pointed out today an accelerated expansion of the uni- 
verse, dominated by an exotic fluid of negative pressure. Furthermore, there are 
evidences suggesting this exotic fluid could be of phantom nature [11]. Hence, 
several classes of black holes (also wormholes, see [12]) have been found which 
have phantom characteristics. An important class of phantom solutions is the 
Gibbons and Rasheed's of the EMD theory [13]. Several others generalizations 
were obtained, such as the higher-dimensional black holes by Gao and Zhang 
[14] and the higher-dimensional black branes by Grojean et al [J '>]. The analysis 
of the algebra produced by a metric with two times in higher dimensions, which 
provides phantom fields in 4D, was developed by Hull [16] and for Sigma models 
by Clement et al [17]. On this work, we will study some solutions coming from 
the EMD theory, which were studied in detail in [18]. 

There are some methods of analysis in black hole thermodynamics theory 
dubbed as geometrical, because they make use of differential geometry to de- 
termine thermodynamic properties such as: points at which black holes become 
extremal or they pass through a phase transition and thermodynamic stability 
of the system. One of the first methods was proposed by Rao [7] , subsequently 
developed by other authors [8]. Later, the work of Weinhold [0] and Ruppeiner 
[10] were frequently used for the study of the black hole thermodynamics. The 
method we will explore in this work is known as Geometrothermodynamics 
(GTD), however, the results obtained with this method will always be com- 
pared with the ones obtained by the usual non-geometric methods. The GTD 
had been shown to be equivalent, and in some situations, even superior in many 
aspects when compared with the usual non-geometric approaches. The GTD 
has been widely used in the literature to study the most diverse classes of black 
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holes [19, 20]. The results obtained by the GTD methods reconciles some incon- 
sistencies between the Weihold and Ruppeiner methods, for example, we can 
mention some cases in which a black hole suffer a phase transition according 
to one method but not by other one. Yet, there is still a good concordance 
when different thermodynamic potentials are chosen, like for example the mass 
and entropy representations, this in virtue of the invariance of the formalism 
by Legendre transformations. Finally, we have mention that the results are 
independent of the particular thermodynamic ensemble considered. 

However, we also have to point out that the GTD method can contain some 
inconsistencies when compared with the more usual analysis done by the specific 
heat. Recently it has shown that for the cases Reissner-Nordstrom-AdS and 
(phantom case) anti-Reissner-Nordstrom-AdS black holes [21], the GTD method 
does not reproduce the results obtained by the specific heat method. We will 
arrive to the same conclusion here in the case of phantom black holes within 
the EMD theory. 

This paper is organized as follows: In Section 2 a summary of the static and 
spherically symmetric solution for the EMD theory in 4D 3 and the derivation 
of thermodynamics variables are presented. In Section 3 the GTD method is 
reviewed in some detail. Section 4 is divided in two parts, the GTD method is 
applied to the classes of asymptotically flat black holes of the EMD theory and in 
particular to anti-Reissner-Nordstrom case in sub-section 4.1, and the analysis 
of the local and global stability is developed in sub-section 4.2. Conclusions and 
perspectives of the work are presented in Section 5. 

2 Phantom black hole solutions and the first 
law of thermodynamics 

In this section we present the class of solutions coming from Einstein- 
Maxwell-Dilaton theory, their relevant parameters, the formulation of the first 
law of thermodynamics, as well as the fundamental ingredients necessaries for 
a detailed analysis of the thermodynamics properties of these solutions. 
We begin by defining the Einstcin-Maxwcll-Dilaton (EMD) action as 

S=f dx A ^[n~2i ll g^V^V v ^ + r ]2 e 2 ^F^F llv \ , (1) 

where the first term corresponds to the usual Einstein-Hilbert action, the second 
one is the kinetic term of the scalar field (dilatonic or phantom) and the third 
one is the coupling term between the scalar and Maxwell fields, with real valued 
coupling constant A. The coupling constant 771 takes the values 771 = +1 (dila- 
ton) or rj\ = —1 (anti-dilaton) . The parameter 772 can be 772 = +1 (Maxwell) or 
7/2 = — 1 (anti-Maxwell). Depending on whether the contribution of the energy 
density is positive or negative, which is determined by 771.2, the kinetic term 

3 From here it will be implied that the black hole solutions discussed here are static and 
spherically symmetric solution coming from the EMD theory in 4D. 
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of the scalar field and the coupling term with the Maxwell field can be normal 
(ordinary) or phantom. 

Now, we use some well established results about the class of black hole 
solutions in EMD theory. According to [18], Eq. (2.19) in Section II, let us 
choose the solutions of u>(u) and J(u) as the function sinh(w), where the horizon 
is non-degenerate and uq > 0, and let us perform a reparametrization of the 
radial coordinate as 

1 / ± = 1-^, (2) 



(r+ - r_) V/-/ r 

with 

r , 1- 

' ^ _ g=p2au 



r ± =± , 2 L n ,„ (r+-r_=2a), (3) 



then we obtain the solution 

dS 2 = f+fldt 2 f+'f^dr 2 -r 2 /^W , (4) 
F = -^drAdt , e" 2Av = , (5) 

where tpo — 0, 7 = A_/A + (for 771 = 1, — 1 < 7 < 1, and rji = —1,7 € 
(-oo,-l) U (l,+oo)), < r_ < r + for rj 2 ^+ > 0, and finally r_ < < r + 
for r)2\+ < 0. This is the exact solution of a spherically symmetric black hole, 
asymptotically flat, electrically charged and static, with internal horizon "r_" 4 
and event horizon which is related to the physical parameters, mass and 

charge of the black hole, trough the relations 

M = + 7 , (6) 



1 + 7 

— — y/mr+r- ■ (7) 

Now, we are interested in the geometrical analysis representing semi-classical 
gravitational effects of the black hole solutions mentioned before. By semi- 
classical we mean quantize the called matter fields, while leaving classical the 
background gravitational field. Therefore we will work with the semi-classic 
thermodynamics of black holes, studied first by Hawking [1], and further devel- 
oped by many other authors [28]. 

There are several techniques to derive the Hawking temperature law. For ex- 
ample we can mention the Bogoliubov coefficients [ ] and the energy- momentum 
tensor methods [4, 28], by the euclidianization of the metric [29], the transmis- 
sion and reflection coefficients [31, 32], the analysis of the anomaly term [38], 
and by the black hole superficial gravity [ ]. Since all these methods have 
been proved to be equivalents [37], then we opt, without loss of generality, to 
calculate the Hawking temperature by the superficial gravity method. 

4 A detailed discussion of the causal structure for the phantom case can be found in [18]. 
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The superficial gravity of a black hole is given by the expression [35] : 



9oa 



(8) 



r — r+ 



where r + is the radius of the event horizon. The relationship between the 
Hawking temperature of the black hole and its superficial gravity is given by 
the formula [1, 34]: 

T= h- (9) 

Therefore, for the the black hole defined by (4) , we have that the superficial 
gravity (8) takes the form 

_ 0+ - rJ) 1 



and the corresponding Hawking temperature (9) would be 

(r+ - r_) 7 



T : 



47IT 



,1+7 



(10) 



(11) 



We define the area of the horizon of the black hole as 



A = 



V922933 dOd(j) 



= 47rr 2 /-~ 7 



= 47rr^+ 7 



(r+ - r_) 



1-7 



Then the entropy of the black hole can be defined as [ ] 
S= ^-A = irr^(r + -r_)^ 7 . 



(12) 



(13) 



On the other hand, working out (5) we obtain the electric potential at the event 
horizon which reads 



A = 



F 10 (r')dr' 



+ 00 



r— r_j- r_)_ 



(14) 



Using eqs. (6), (7) and (13), we can write the differential forms of the mass, 
charge and entropy: 



dM = \ (dr + +7dr_) , 

d i = m (rit) ( r - dr + + r + dr -) 



' 1+7 

dS= 



^ { [2 - (1 + 7)^] dr+ - (1 - 7 )dr_ } 



(15) 



these relations, together with (11) and (14), satisfy the first law of black hole 
thermodynamics [2] 

dM = TdS + rj 2 A dq . (16) 

Note that when 772 = —1, the first law is generalized for the Einstein-anti- 
Maxwell-Dilaton sector, where the second term, which is related to the work, 
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suffers a change of signal as a consequence of a negative energy contribution 
to the system. The exact formula of Eq. (16), originally known as the Smarr 
formula [39], can be integrated resulting in 



M = 2TS + m A q . 



(17) 



In this case we also have a change of signal of the work term for the phantom 
case. 

We can write the horizons in terms of the mass and charge parameters of 
the black hole as follows: 



It is important to note that the possible extremal case, i.e., r + — r_, have 
to be analyzed carefully, because according to the causal structure described 
in [18], for the ordinary EMD theory r_ is a curvature singularity, and for the 
phantom case in which r)2\+ < (r_ < 0) is a true singularity . Hence, in the 
last case we have a causal structure identical to the Schwarzschild solution with 
just one horizon. In some phantom cases where 7 takes discrete (integer) values, 
we have two horizons that can be crossed where the radial coordinate r for the 
internal horizon is no longer r_, but r = 0, with r_ being a singularity where 
the geodesies ends. When this is so, it is not simple to determine if the extremal 
limit exist, because r + — > represents a regime which has not yet been analyzed 
in the literature. The analysis by the Geometrothermodynamics method could 
provide a new insight to understand this pathological solution, but we will see 
that, just as the usual method using the specific heat, some subtleties are still 
unavoidable. Maybe this subtleties can be understood only within the context 
of more fundamental quantum analysis. 

When the normal EMD theory is considered, the extremal limit still provides 
a structure quite similar to what is called a naked singularity; classified as a 
light-like naked singularity, which can be reached only after an infinity lapse of 
time. A semi-classical analysis of such a structure, i.e., a non-asymptotically 
flat black hole, has been studied in [31], but the analysis does not contain a 
conclusive result in relation to the thermodynamic properties. So once again it 
seems like some fundamental quantum theory of gravity should be considered 
at this point to depth into this question. 

We have briefly discussed here the causal structure and the extremal limit of 
this class of solutions, now we are interested in their thermodynamic properties, 
which we will see present subtleties precisely for the particular cases considered 
in this work. 

Taking the eqs. (13), (14), (18) and (19), we can rewrite the temperature, 
entropy and the electric potential as functions of the mass and electric charge 




(18) 



(19) 
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as follows: 



T(M,q) = 




1+7 



(20) 



(7-l)M+(7 + l)WM2- 



S(M,q) 



(1 + 7) 



-1-7 



(21) 



( 7 -l)M+(7 + l)WM 2 - 



2?7279 2 
(1 + 7) 



1-7 



A (M, g ) 



79 V V (1 + 7) 



(22) 
(23) 



Later on we will use the entropy as thermodynamic potential. If we want to 
consider the mass, which is equivalent to the energy, as thermodynamic potential 
we need to invert (22) in order to write the mass in terms of the entropy and 
the electric charge. This is by no means evident from the general case (4). 
Thus if we want to proceed with the analysis in terms of at least two different 
thermodynamic potential, i.e., entropy and mass, it is convenient to specialize 
to the case A = 0, 7 = 1 (with r\\ = 772 = 1), which is known as the Reissner- 
Nordstrom, or (—171 = T}2 = — 1) which would be the anti-Reissner-Nordstrom. 
For these cases, the mass, electric charge, entropy and electric potential are 
given by 



M=|(r + + r_), q = ^/rj 2 r + r. 



T = 



S = 



A = ^ 



(24) 



and satisfy the first law of black hole thermodynamics (16) and the Smarr 
formula (17). 

Finally, from (24) we have 



' + 



M + \/M 2 - i m 2 , r_ = M — ^[M 2 



mi 



(25) 



which gives us the temperature, entropy and electric potential which explicitly 
read 



2 ' 



(26) 



2ir (M + y/M 2 - 772 g 2 ) 
S = 7T (M + v/M 2 -r/ 2(? 2 ) 2 , A Q = - (M - v/M 2 - ry 2(? 2 ) /q . (27) 
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We have introduced all the necessary prerequisites to start with the study of 
the thermodynamic properties by means of the framework of Geometrothermo- 
dynamics. In the next section we will introduce this method and further apply 
it to our particular cases of interest. 



3 The Geometrothermodynamics method 

The Geometrothermodynamics (GTD) make use of differential geometry as 
a tool to represent the thermodynamics of physical systems. Let us consider the 
(2n + l)-dimensional space T, which coordinates are represented by the ther- 
modynamic potential $, the extensive variable E a and the intensive variables 
I a , where a — 1, ...,n. If the space T has a non degenerate metric Gab(Z c ), 
where Z c = £ a , I a }, and the so called Gibbs 1-form 9 = rf$ - 5 ab I a dE b , 
with S a b the delta Kronecker; then the structure (T, O, G) is said to be a contact 
riemannian manifold if O A (dQ) n ^ is satisfied [30]. The space T is known as 
the thermodynamic phase space. 

We can define a n-dimensional subspace E C T, with extensive coordinates 
E a , by the map tp : E -> T, with $ = <$>(E a ), such that 

I dW - 6 *b± ■ 

We call the space E the thermodynamic space of the equilibrium states; 
the first equation (28) is called "First law of thermodynamics" , and the second 
relation would be referred as the "condition of thermodynamic equilibrium" . 
We impose as a necessary condition the "Second law of thermodynamics" : 

± dE^E^^ ™ 

where the signal (±) depends on the chosen thermodynamic potential, for ex- 
ample, in the case of the mass we have (+), and for the case of entropy we 
have (— ); this is known as convexity of the thermodynamic potential condition. 
The thermodynamic potential is defined such that it satisfy the homogeneity 
condition Q(aE a ) = a^^?(E a ). By differentiation with respect of a, using the 
second equation of (28) and taking a = 1, we have that 

P$(E a ) = 5 ab I a E b . (30) 

By differentiation of this last equation, and using the first equation of (28), 
we have 

(1 - P)8 ab I a dE b + S ab E a dI b = . (31) 

For (3 = 1 we get the Euler identity (30) and the Gibbs-Duhem formula (31). 
The pullback tp* : T*(T)(g>T*(T) ->■ T*(E)®T*(E) 5 , induce a metric on E, such 
that ip*(G) = g. 



5 Where T*(T) and T*(E) represent the tangent spaces of T and E, respectively. 
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Hernando Quevedo has developed and improved a possible metric G for the 
GTD, which for the case of black holes it can be written as [(>]: 

dL 2 = G AB dZ A dZ B = 6 2 + (5 ab E a I b ) (r] ab dE a dI b ) , (32) 

where rj ab = {±1, L, •••! !}• The case r) ab = {— 1, 1, 1} holds for second order 
phase transition. The metric on E induced by the pullback ip* is: 

87 A f)Z B 

dl 2 = g ab dE a dE b = -^^-rG AB dE a dE b 
y dE a 8E b 

= { E '§;){^ s, 'iS^) dE ° dE> - < 33 > 

Since the thermodynamics system has to be independent of the particular 
choice of thermodynamic potentials, and also invariant under Legendre trans- 
formation, we have that the metrics (32) and (33) should be invariant under 
Legendre transformations of the form 

$ = $ - S ab E a 7 h , E a = -I a , I b = E b . (34) 

The n-dimensional space E, with metric g ab , contains the information about 
the thermodynamic interaction, phase transitions and fluctuations or stability 
of the thermodynamic system. With the help of metric (33), we can derive 
the scalar curvature R, which gives us information about two aspects of the 
theory: when there are thermodynamic interaction and when there are phase 
transitions, also it tell us at which points of the thermodynamic equilibrium 
space these transitions take place. 

As we mentioned in the introduction, this method presents some important 
results. The first one is that the use of the Weinhold and Ruppeiner metrics 
provide results which are in contradiction to each other, even in contradiction 
to themselves when, for example, different thermodynamic potentials are used 
to describe the same system. We can cite the case of the Reissner-Nordstrom 
(RN) black hole for which the use of the Ruppeiner metric in the entropy repre- 
sentation provides a flat space E, hence without phase transition [23]. However, 
when the representation is given by the internal energy of the black hole, the 
same method points out a non-zero scalar curvature with a singularity, i.e., 
with a phase transition in the curved space E [22]. This contradictory result 
is resolved in the GTD by the metric (32), which is invariant under Legendre 
transforms (34), reconciling the results of RN black hole thermodynamics re- 
gardless the choice of thermodynamic potential [24]. The second result is the 
agreement with the usual analysis of the thermodynamic system by means of 
the specific heat of the black hole [20]. A third result is that the description 
of the thermodynamic system depends, in many cases, of the ensemble's choice 
[26, 25], because this choice leads to different specific heats. This problem is 
solved by the description of GTD, so it results in a consistent description of the 
thermodynamic system described by different ensembles [27]. 

We saw that the metric of the thermodynamic equilibrium space E, can 
be obtained by the pullback of the metric defined on the contact riemannian 
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space. By the definition of the line element (33) in E space, we can define the 
distribution of probabilities to get a physical state with extensive variable E a 
within the interval E a + dE a : 



P(^) = ^Pexp 



9ab dE a dE b 



where P{E a ) satisfy 



/n 
J[ dE a P{E a ) = 1 



(35) 



(36) 



It can be shown that by taking the derivative of (35) with respect of V , 
where V is the volume of the system, we obtain the expression for the second 
fluctuations (in the thermodynamic limit V — > oo) [40] 



(AE a AE b 



-9 



ah 



(37) 



where AE a = E a - £ ( a 0) , and g ab is the inverse of g a b- A more realistic analysis 
requires the constants to be adjusted. If the fluctuations are small and real 
valued, then we say the system is stable. 

Another criteria used to determine the stability is trough the following geo- 
metric objects 



P? 



g n > , pf ] = g 22 > , 



Pi 



(n) 



9nn > , 



(38) 



9n 


912 


912 


922 


9ii 


912 


912 


922 


913 


923 



922 923 

923 933 



>0, 



> , p n = det [g a b] > 



(39) 
(40) 



The positive (negative) signal of p n , which would depend on the choice of the 
thermodynamic potential, determines the local (un) stability of the thermody- 
namic system; more specifically if we have 



Pi > 



1. 



(41) 



then we can affirm the system is globally stable. 

Since we are not considering rotating black holes, another way to determine 
the global (un) stability is through the Hclmholtz free energy. In terms of black 
hole thermodynamics variables, the Hclmholtz free energy is no other thing than 
the Legendre transformation of mass (energy) M(S,q): 



F(T, q) = M(S, q) - TS 



(42) 
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When we have 



F(E a ) < , VE a (E a G I(E a )) , (43) 

where I(E a ) is an interval and E a are the extensive variables, then the ther- 
modynamic system is said to be globally stable. In the case of usual thermody- 
namics the Helmholtz free energy is given by F(T, V) = U — TS. 
We can also define the Gibbs potential as 

G(T,A )=M(S,q)-TS-ri 2 A q, (44) 

in this case the global stability is determined by 

G(E a ) < , VE a (E a G I{E a )) . (45) 

We introduced the signal 772 in (44) to compensate the contribution of the work 
term. Later we will make use of the Gibbs potential to determine the global sta- 
bility of the thermodynamic system. For a system shown stable in this ensemble, 
we must have 

d 2 G d 2 G d 2 G 

dfi'dA 27 OTdAo ~ ' ( ' 

In the next section we will apply these methods to determine the thermody- 
namic properties of the black hole solutions coming from the EMD theory. 



4 Thermodynamics of phantom black holes 

4.1 Application of the Geometrothermodynamics method 
4.1.1 Einstein-(anti)Maxwell-Dilaton solutions 

To begin with, let us define the thermodynamic variables of the system. 
For the Einstein-Maxwell-Dilaton black holes, we will have always a solution 
with two physical parameters, the mass M and the electric charge q. Other 
variables (parameters) such as entropy S, temperature T and electric potential 
Aq, can be defined as implicit functions of the parameters mentioned before. 
The contact riemannian manifold T is, in this case, 5-dimensional and the space 
E of thermodynamic equilibrium states is a 2-dimensional manifold. 

The thermodynamic description is the entropy representation, S(M, q), which 
is identified as the thermodynamic potential $, according to was defined in the 
previous section. The extensive variables are the mass M and the electric charge 
q, which are represented by the coordinates E a . The intensive variables are the 
temperature T and the electric potential A , represented by the coordinates I a . 

We have then a coordinate system for the thermodynamic phase space T as 
being Z A = {S(M, q),M, q, T, A }, together with the Gibbs 1-form given by 6 

Qs = dS-~dM+ r]2 ^dq, (47) 

6 This expression comes from the first law of thermodynamics (16), which was inverted in 
order to isolate the entropy. 
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such that ip* (&s) — is satisfied, which is no other thing than first law of black 
hole thermodynamics: dM = TdS + ^A^dq (772 = ±1). 

For a second-order phase transition, the line element (32) of the space T is 



dL 2 



dS dM 

T 



V2A0 
T 



dq 



-d\ -\dM + d 



T 



M 772A0 
T 

dq 



T 



q x 



(48) 



The first and second laws, and the equations of the equilibrium state as are 
given by 



d§ = S ab I a dE' 

d 2 S d 2 S d 2 S 



dS = ±dM-^dq 
T T 



dM 2 ' dMdq ' dq 2 



sS 



9$ b dS _ 1 dS 



T 



(49) 
(50) 
(51) 



Now we have to specify a solution, where we consider first the general case 
(4). The line element (33) of the equilibrium space, taking into account (22), 
would be 



„., , , , dS dS 

,lh = {M dM +q dq 



d 2 S Jnr2 d 2 S l2 
dM 2dM + W dq 



guudM 2 + g qq dq 2 , 

(m + ^m 2 - 2 ^) 



(52) 
(53) 



1+27 



9mm = 167r 2 7 27 -- 



M 2 2^ 



x M 2 — rj 2 q 2 +Ma M 2 - 



(1 + 7) 



( 7 -l)M+( 7 + l)JAf 2 - 



(1+7) 



-1-27 



(1+ 7 )M-(1 + 3 7 WM2-^) 



(54) 



9qq 



-32ry 2 7T 7 



2^,1+27 



M 



x M 2 



V M2 


2?7279 2 \ 
(1+7) ) 




2jj279 2 
(1+7) 


M 2 - 


2r) 2 jq 2 \ 
(1 + 7) J 



-1+27 



(55) 
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-1-27 



(56) 



We notice that the phantom contribution can switch the signature of the metric 
of the space E. To calculate the scalar curvature associated with the metric 
(53), we can make use of a particular formula valid for 2-dimensional spaces 



R(M,q) = - 



1 



V\ det l9]\ 
det[H s ] 

2{det[g]f 



(dqQMM — dklQMq I q I t>M9qq ~ 9 q gMq 
7HH J w 



V V\ det l9] 



9MM 9Mq 9qq 

9m9MM 9M9Mq 9M9qq 

d q 9MM d q gM q 



(57) 
(58) 



Replacing (53) into (57), we have 

R(M,q) 

where 



N(M,q) 
D(M,q) 



(59) 



N(M, q) = S 2 MM {FFqSqqq + 2S qq {F 2 q - FFqq)) 

+ FSqq { — SqqFMSMMM + F {S\ [Mq — SuMMSqqq)) 

+ Smm{ — 2SqqF M + FSqq { — FqSuMq + FklSqqM + 2S qq FMM) 



F 2 (S 



with di...djS 



where 



MMqSqqq 



Si...!, n,...<),F 



n2 

°qqM 



2S qq (SmMMM 

Fi,„j, and 



S, 



qqM M 



)) 



D(M,q) = 2F z S 2 rM S 2 



'MM^qq 1 



F = MS 



M 



qSq , 



(60) 

(61) 
(62) 



and S as given before by (22). 

With the help of a mathematical software, we find the zeros of the scalar 
curvature R(M, q) (see (59)) are the mass values Mi — ±g^/2?727/(l + 7), 
M 2 = ±9^/^2(1 + 1)/2- We also have another zero for the scalar curvature 
for the particular case 772 = —1 at M 3 = g-y/(l + 7)/2(l + 27). This tell us that 
in general the scalar curvature is non zero - which implies the existence of ther- 
modynamic interaction for this class of black holes [41] - and is zero only when 
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the mass takes the values Mi, M 2 or M3. These values characterize when the 
the black hole becomes extremal. Let us note that only when we set r + = r_ 
in (18)-(19), we obtain the value M 2 for the mass. Probably this is related to 
the fact we pointed out before, that in the normal and phantom cases, r__ is a 
true singularity of the curvature, therefore, the extremal limit in this regime is 
not viable. The two values Mi and M 2 , for 7y 2 = 7 = 1 (RN), are Mi )2 = ±q, in 
total agreement with [4] . We also want to note that there exist an extremal limit 
at M 3 , when 772 = —1. The point at which the scalar R(M,q) diverge is given 
by the mass M 4 = ±q{l + 3"f)/ y/2r) 2 (l + 7)(1 + 27). Then, when ry 2 = 7 = 1, 
M4 = ±2q/y/3, also in agreement with the previous work of Davies [4]. How- 
ever, for the anti-RN solution M4 is not real, which means that there is no phase 
transition for the phantom when work with the entropy representation. 

To corroborate our analysis, let us calculate the specific using the well know 
formula 



C a 



(dM 
\~dT 



(dM 
\dS 



d 2 M 
~d~S 2 ~ 



( dS_\ 2 I ( d 2 S\ 



(63) 



Doing this we obtain: 



C„ 




4tt7 7 Jm 2 - 



(7-l)M+(l + 7)^M2-to 



+7) 



(1 + 7 )M-(1+3 7 )^M2-^ 



1+7 



(64) 



The zeros of (64) give us the information about the points at which the black 
hole becomes extremal; these are precisely Mi and M 2 founded by the previous 
analysis, thus corroborating our analysis. In the same way, we get that the 
phase transition point at which C q diverges is M4. Remarkably is the fact that 
M3 is not a zero of the specific heat, hence, it should be some spurious zero that 
results from a failure of applicability of the analysis for this class of pathological 
solutions. We will see below that not only a new (non physical) extremal case 
is revealed by the GTD method, but also a new critical point is revealed when 
we choose the mass of the black hole as the thermodynamic potential. 



4.1.2 (anti)Reissner-Nordstrom solution 

To study in detail the thermodynamics of the EMD system, let us consider 
the particular case 7 = 1, this corresponds to the (anti)Reissner-Nordstrom 
((anti)RN) solution. For this case, we have that the entropy is 

S = 7T (m + ^M 2 -V2q 2 ) 2 , (65) 
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that, when replaced the into (53), furnish the following: 



dl 2 



M 2 - m q 2 + My/ M 2 - m q 2 
(M 2 - m q 2 ) 2 



x{(M+ y/M 2 -mq 2 ) (-M + 2y/M 2 - r] 2 q 2 ) dM 2 
+ 772 (M + y/M 2 - T72-7 2 ) [(M 2 - mq 2 f /2 + M 3 ] dq 2 } . (66) 
The scalar curvature derived from this metric is given by (59), where 



N(M, q) = -y/M 2 - m q 2 [M + y/M 2 - r] 2 q 2 ) 



1024M 14 - 3840?72(7 2 M 12 



5504g 4 M 10 - 3408t7 2 <? 6 M 8 + 460g 8 Af 6 + 377r7 2 g 10 M 4 
122g 12 Af 2 + 5?72q 14 + \/M 2 - r;2g 2 (l024m 13 
3328?72g 2 A / / 11 + 3968<? 4 A/ 9 - 1776rj 2 q 6 M 7 - W0q 8 M 5 



251r;2Q 10 M 3 - 35g 12 M 



(67) 



D(M,q) = 4tt 5 



64M 10 - 208r] 2 q 2 M s + 248g 4 M b - l37r) 2 q b M 



.2 fl/ r8 



,6 A/f4 



34g 8 M 2 - 2T] 2 q W + V ' M2 - mq 2 X f 64M 9 - 176?72<Z 2 Af 7 



168g 4 7\f 5 - 7lr ]2 q M' s + llg s M 



6 ]\j3 



(68) 



The zeros of the numerator, N(M,q), in (67), are M5 = ±q^/fj2 and Mq — 
±iqy/r] 2 /3, which depend on the choice of 772- The RN case, when 772 = 1, has 
the zero M5 = ±q (where we did r + = r_ ) that corresponds to the extremal 
RN black hole. The anti-RN case, 772 = — 1, has the zero Mq = ±q/>/3, but this 
result reveals a weakness of this method, because it points out the presence of 
an extremal anti-RN black hole, which we know does not exist. The analysis of 
the causal structure performed in [18] shows that the anti-RN black hole has a 
causal structure identical to the Schwarzschild black hole, hence, there is no exist 
extremal limit for this case. We don't have any explanation for such irregularity 
of the thermodynamic system when described by the GTD method. We believe 
that the extremal limit for non trivial black holes, just like the anti-RN case, 
reveals pathologies which are not well described by the GTD. 

The zero of the denominator D(M, q), in (68), is given by M-? = ±2qy/ri 2 /3, 
that is real only for r/ 2 = 1. It occurs then that the RN black hole has a second 
order phase transition point at M 7 , in good correspondence with Davies [4]. On 
the other hand, the anti-RN black hole, 772 = — 1, has not any phase transition 
point; consequently, there is not extremal analogue, nor phase transition for the 
anti-RN case. 

As before, order to check our results, let as can calculate the specific heat of 
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these model. The specific heat (63), calculated in combination with (65), is 



_(m + Vm 2 - mq 2 ) 

C q = 2tt^AP- r, 2 q 2 ) J . (69) 

M - 2y/M 2 - r\ 2 q 2 



From (69) we have that the black hole is extremal (C q = 0) when M = M$ = 
±qy/f]2, and it has a point of phase transition {C q — > oo) at M = M 7 — 
±2qy/rj 2 /3 . These results are in total agreement with the previous results we 
found by the GTD method, but we want to stress that before we obtained an ex- 
tra extremal anti-RN black hole solution. This reinforce our previous statement 
about the weakness of the GTD method for some pathological situations. 

As was mentioned before, we will consider now the mass representation anal- 
ysis. We start by inverting (65) to write the mass in terms of the entropy and 
electric charge, this is 

M{S t q)=S±^£. (70) 
The Gibbs 1-form reads in this case 

9 M = dM - TdS - r) 2 A dq . (71) 
Then, using (32) we obtain 

dl? = {dM - TdS - r) 2 A a dq) 2 + (TS - r] 2 A q) [-dTdS + d (-r/ 2 A ) dq] . (72) 
The pullback tp* induce a metric on the space E, 

„ a ( Q dM dM\( d 2 M 2 d 2 M 2 \ 

dl =iS ^S +q ^q-){'^ dS + ^ dq ) (73) 



(S + 3r, 2 Trq 2 ) 
4S 



(S-3 V2 nq 2 ) 
8irS 2 



dS 2 + r, 2 dq 2 



(74) 



Finally, with this metric we can compute the scalar of curvature, 



rtf a \ 2^ 2 q 2 S 2 {S- m 2 ) 

R{S,q) = r] 2 — — ^3 . (75 

We can see that when the thermodynamic potential is the mass of the RN 
black hole, the interpretation by this method of the extremal limit (r + = r_) 
is clear: the zeros of the numerator of R(S, q) in (75) exist only for the values 
Si = irq 2 , which clearly represent the extremal RN black hole. When the 
anti-RN (phantom) case, rj 2 = — 1 , is considered, we evidence that there is no 
extremal limit, as it should be. This last conclusion solve in some sense the 
unphysical prediction we get in the entropy representation, by changing to the 
mass representation. But in contrast a new problem arise, we get a new point 
of phase transition for the anti-RN case at 5*2 = 3rj 2 irq 2 , with r + = 3ry2?'-. This 
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indicates a breach of the invariance of the theory for describe the system regard- 
less of choice of thermodynamic potential. Once again this result confirms our 
claim about the failure of the GTD method when applied to pathological solu- 
tions, which despite having a simple causal structure, apparently they present 
non trivial solutions. For the RN solution we have a extremal limit at Si and a 
second order phase transition point at S2, in agreement with [4]. 

Finally, let us compute the specific heat. Using (63) together with (70), we 

From this formula the correct interpretation can be read: for the RN black hole 
we have an extremal limit when S = Si = irq 2 , and a phase transition point at 
S = S2 = 3irq 2 . On the other hand, for the anti-RN case, we conclude there is 
no extremal limit, nor point of phase transition. 



4.2 Local and global stability- 
It is usual to study local stability of a thermodynamic system by means of 
the specific heat. Alternatively, within the GTD, we can study the metric com- 
ponents of the thermodynamic equilibrium space E or even more, the Hessians 
of the entropy and mass. To determine the global stability it can be achieved by 
the analysis of all the components of the metric as well as their corresponding 
determinants; but also by means of the Helmholtz free energy or by the Gibbs 
potential. Here we will continue with the study of local and global stability of 
EMD black holes solutions. 

Let us start calculating the Hessian of the entropy for the general case (22), 
which is defined as 

/ d 2 S d 2 S \ 

Hs=l %_ d ^ q ■ (77) 

\ dMdq dq 2 ) 

Using (22) and considering 772 = 1 (including 7 = 1) we get local instability in 
view of Smm(M, q) is always positive, and S qq {M, q) and Suq{M, q) are positive 
as long as q > 0. In the same manner, when 772 = —1, Smm(M, q) and S qq (M, q) 
are always positive and SM q (M, q) is positive provided that q > 0; therefore we 
have that the solutions are locally unstable. 

Analogously, we can calculate the Hessian matrix of the mass for the partic- 
ular choice 7 = 1. Using (70) we found that 

( d 2 M d 2 M \ 
a 2 M 9 a 2 M ( 78 ) 
dSdq dq 2 ) 

C 3r]2Trq 2 — S y / 7ri)2g 
8v=S 5 /2 2S 3 / 2 
2S3/2 Sl/2 

The Hessian of the mass also leads us to local instabilities. More precisely, when 
r)2 = 1, then M$s and M$ q can take positive or negative values. Similarly, when 
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?72 = —1, Ms q can takes positive or negative value, whereas M$s and M qq are 
always negative. 

Once again, let us consider the analysis by the study of the specific heat 
(64), which in terms of r + and r_ can be written as 



(r+ - 7 r_)(r_ ( 



\l-7 



[r+ - (1 + 2 7 )r_ 



(80) 



From this expression, we can establish, for example, local stability (C q > 0) 
for the event horizon interval 7r_ < r + < (1 + 2 7 )r_, < 7 < 1 (EMD with 
rji = 1). The specific heat, with 7 = 772 = 1 or 7 = —772 = 1, corresponds to 
the RN or anti-RN cases, respectively. The anti-RN case with 7 = 1 is locally 
unstable because r_ < implies C 9 < 0. The only case where the phantom 
solutions are locally stable will be when 7r_ < r + < (1 + 27)r_ for 7 < — 1. 

As we mentioned before, the components of the metric of the thermodynamic 
equilibrium spac, E, give us also information about local stability. According to 
(52), the components of Qmm and g qq , written in terms of r + and r_, read 



9mm 



9qq 



-167rV (1+7) 



o 2 I+27 



[r+ - (l + 2 7 )r_] 
(r + — 7r_)(r + — r_) 2 T ' 
[r^ + (7-l)r_+2 7 (7+i)r 2 ] 
( r + ~ 7 r -)( r + ~ r-) 27 



(81) 
(82) 



These components always have opposite signals for each of the two cases, normal 
and phantom, including the case 7 = 1. Therefore, by simple inspection of the 
components of the metric on E, we assert that the system is locally unstable; 
additionally, it tells us that there is global instability, just as was mentioned in 
Section 3. 

For the next step, let us study the global stability of the class of solutions 
of interested. For the general case we can write 



M = 



q 


H 


2r?27^ol 


2A 


(1 + 7). 


1 


1 - 




4ttT 


(1+7) J 



4A 



1 - 



2%^ 1 



(1 + 7) 



(83) 
(84) 



The analysis start with the grand canonical ensemble. Using (83) and (84), we 
calculate the Gibbs potential (44): 



G(T,A ) = 



1 



16ttT 



1 - 



(1 + 7) 



1+7 



(85) 



We know that T > 0, so we have to analyze the term between brackets. 
Consider first the normal case with 772 = 1; when (1 +7) is an odd integer, then 
the system is globally stable only for 7 > and A G (—00, — y(i + j)/2) U 
(iy/ (1 + 7)/2, +00). The critical electric potential values are A ( c \ = (1 + 7) 
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(y/2 ) (by setting r + = in (84)). No let us consider the phantom solutions 
with r]2 = — 1, in this case the system is globally unstable, except for 7 < — 1 and 
Aq 6 (— 00, — yjl + 7|/2) U 1 + 7|/2, +00). A straightforward calculation 
of the derivatives in (46) shows that 



d 2 G 
dT 2 

d 2 G 
dAl 
d 2 G 



1 



8ttT 3 
(1 



(1 + 7) 
2V2A 2 ) 



1+7 



4ttT 



2^ 
(1 + 7) 



7-1 



1 - 



2^ 
(1 + 7) 



(86) 
(87) 
(88) 



From here we can see that (88) always break the stability criteria, because 
Aq can be positive or negative (remember that q is real). Then it follows that 
the system is always locally unstable. 

The corresponding analysis by the canonical ensemble formalism is slightly 
different. We start writing the Helmholtz free energy (42): 



F = 



1 



16vrT 



1 - 



(1+7) 



1 + 2 m A 2 



l + 2 7 
1 + 7 



(89) 



The analysis of this formula is the following: the system is globally stable for 
two cases. The first ones when A e (— 00, — + 7)/2) U + 7 )/2, +00) 
for 7 > odd integer and i]2 = 1. The second case is when Ao 6 (—00, 
-V( 1 + 7)/2(l + 2 7 )) U + 7)/2(l + 2 7 ), +00) for 7 even integer and tj 2 

I or A € (-00, -^(1 + 7)/2) U (^(1 +7)/2, +00) for 7 e 3?. When we use 



the specific heat formula (80), we found second order phase transition points 
given by A 0(1) = ±y/rj 2 {l 4 
for 7 > 1) . The point A ( 3 ) 



- 7 )/2(l + 2 7 ) and A 0(2) 



7)/2 ( even 

7)7(27), as well as A (2) with 7 < 1, 



represent the extremal case where r + — > (see [18] for details). 

Finally, we can calculate the minimum temperature for this system using 
the formula dT~ l /dr + = [42], so we have: 



T 



1 



7 



(90) 



47rr + y 1 + 7, 

From the specific heat, Eq. (80), we identify a phase transition point at r + = 
3r_, which gives us the critical temperature 

T, = — ! — I , (91) 



47rr + \ 1 + 27 
combining this with (90) we get the relation T c 



T [2(l + 7)/(l + 2 7 )]- 



5 Conclusions 

The thermodynamic properties of the class of solutions known as phantom black 
holes has not been studied in detail yet. The physical stability of these solutions 
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can be determined also by studying their thermodynamic properties. The main 
objective of the present work was to fill this gap by establishing a detailed 
analysis of this kind of solutions. 

The zeroth, second and third laws of thermodynamics remain unmodified by 
these solutions. However, the first law had to be generalized to take into account 
the contribution of the work done on (or by) the system with the "wrong" sign 
(when compared to the usual case) . The differential form of the generalized first 
law was written in (16) whereas the exact expression was written in (17). 

The use of Geometrothermodynamics as a tool of analysis has proven to be, 
once more, equivalent to the most usual methods, but with an important ex- 
ception that was the phantom case with 772 = — 1. More specifically, we saw for 
the phantom solutions that when we choose the entropy as the thermodynamic 
potential, the method introduce a new value of the mass parameter which could 
be interpreted as a (non physical) extremal black hole limit. If we choose the the 
mass as the thermodynamic potential, then we get a new critical point for the 
system. The new mass parameter and the new critical point have to be consid- 
ered as spurious zeros of the numerator and demininador of the scalar curvature 
of the space E. This can be understood by the fact that EMD solutions switch 
to the phantom sector (i.e. with 772 = —1) by the symmetry transformation 
q 2 — > — q 2 , as was shown in [ ]. In the case of the entropy representation, with 
the choice 7 = 1, that kind of symmetry leads to the appearance of a new real 
valued zero of the scalar curvature which originally was pure imaginary. In the 
case of the mass representation, the same symmetry transformation entails the 
appearance of a new critical point, i.e., a divergence of the scalar curvature. 
When we carried out the study of the thermodynamic properties by means of 
specific heat, we have not found new spurious critical points nor new divergence 
points of the scalar curvature, just as was expected. This result revealed the 
fragility of the GTD method when applied to pathological solutions. 

In regard to the stability analysis, we found that the only possible local 
stability would correspond to the following cases: 7r_ < r + < (1 + 2-y)r-, for 
Vi — V2 — 1 and < 7 < 1 or for 771 = rj 2 — — 1 and 7 < — 1. All the other 
solutions, normal or phantom, has been shown to be locally unstable. The 
global stability can be established for some particular situations that restrict 
the values the electric potential Aq and the parameter 7 can take. 

The perspectives of the present work is to study in detail the subtleties 
that arise from the phantom solutions, with the hope to strength some weak 
aspects of the promising novel method which proved to be the Geometrother- 
modynamics. But also we expect that will be evidenced some physical limits 
to the use of this method. This has already been shown to be true in the cases 
of Reissner-Nordstrom-AdS (RN-AdS) and anti-RN-AdS black holes [ ]. We 
propose elucidate with more details this issue in a forthcoming work. 

Finally, with respect to the global stability of one portion of the class of 
phantom solutions, we expect this to be an indication for the stability of the 
space-time of these class of solutions. That shall also be a topic to be studied 
in a subsequent work. 



20 



Acknowledgement: We are grateful to Gabriela Conde Saavedra for the 
help in the elaboration of the manuscript. M. E. Rodrigues thanks UFES 
for the hospitality during the development of this work and Z. Oporto thanks 
CLAF/CNPq for financial support. 

References 

[1] S. Hawking, Commun. Math. Phys. 43, 199 (1975). 

[2] J. M. Bardeen, B. Carter and S. W. Hawking, Commun. Math. Phys. 31, 
161 (1973). 

[3] C. O. Lousto, Nucl.Phys. B 410: 155-172 (1993); Erratum-ibid. 449: 433 

(1995) . 

[4] P. C. W. Davies, Proc.Roy.Soc.Lond. A 353: 499-521 (1977). 

[5] Jorma Louko and Stephen N. Winters-Hilt, Phys. Rev. D 54:2647-2663 

(1996) . 

[6] H. Quevedo, A. Sanchez, S. Taj and A. Vazquez, Gen. Rel. Grav. 43:1153- 
1165 (2011). 

[7] C. R. Rao, Bull. Calcutta Math. Soc. 37, 81 (1945). 

[8] S. Amari, Differential-Geometrical methods in Statistics (Springer- Verlag, 
Berlin, 1985). 

[9] F. Weinhold, J. Chcm. Phys. 63, 2479, 2484, 2488, 2496 (1975); 65, 558 
(1976). 

[10] G. Ruppeiner, Phys. Rev. A 20, 1608 (1979); Rev. Mod. Phys. 67, 605 
(1995); 68, 313 (1996). 

[11] S. Hanncstad, Int. J. Mod. Phys. A 21, 1938 (2006); J. Dunkley et al., 
Astrophys. J. Suppl. Ser. 180, 306 (2009). 

[12] K.A. Bronnikov, M.S. Chcrnakova, J.C. Fabris, N. Pinto-Ncto and M.E. 
Rodrigues, Int.J.Mod.Phys.D 17:25-42 (2008). 

[13] GW. Gibbons and D. A. Rashecd, Nucl. Phys. B 476, 515 (1996). 

[14] C. J. Gao and S. N. Zhang, arXiv:hep-th/0604114. 

[15] C. Grojean, F. Quevedo, G. Tasinato and I. Zavala, J. High Energy Phys. 
08 (2001) 005. 

[16] C. M. Hull, JHEP 9807 (1998) 021. 

[17] Mustapha Azreg-Ainou, Grard Clment, Julio C. Fabris and Manuel E. Ro- 
drigues, Phys.Rev.D 83:124001 (2011); arXiv:1102.4093vl[hep-th]. 



21 



[18] Gerard Clement, Julio C. Fabris and Manuel E. Rodrigues, Phys. Rev. D 
79, 064021 (2009); arXiv:hep-th/09014543. 

[19] Danny Birmingham and Susan Mokhtari, Phys.Lett.B 697:80-84 (2011), 
arXiv:1011.6654vl [hep-th]; Hernando Quevedo and Maria N. Quevedo, 
arXiv:1111.5056vl [math-ph]; Hernando Quevedo and Alberto Sanchez, 
Phys.Rev.D 79:024012 (2009), arXiv:0811.2524vl [gr-qc]; Alexis Larranaga 
and Alejandro Cardenas, arXiv:1108.2205vl [gr-qc]; Alexis Larranaga 
and Sindy Mojica, Brazilian Journal of Physics 41: 154-158 (2011), 
arXiv:1012.2070vl [gr-qc]; Hernando Quevedo and Alberto Sanchez, Rev. 
D 79, 087504 (2009), arXiv:0902.4488v2 [gr-qc]; Lszl rpd Gergely, 
Narit Pidokrajt and Sergei Winitzki, Eur.Phys.J.C 71:1569 (2011), 
arXiv:0811.1548v3 [gr-qc]; Hernando Quevedo, J.Math.Phys. 48 (2007) 
013506, arXiv:physics/0604164v2 [physics. chcm-ph]; Hernando Quevedo, 
Alberto Sanchez and Alejandro Vazquez, arXiv:0811.0222vl [math-ph]; 
Hernando Quevedo and Alejandro Vazquez, AIP Conf.Proc. 977:165- 
172 (2008), arXiv:0712.0868vl [math-ph]; M. Akbar, H. Quevedo, K. 
Saifullah, A. Sanchez and S. Taj, Phys.Rev.D 83: 084031 (2011), 
arXiv:1101.2722vl [gr-qc]; H. Quevedo, A. Sanchez and A. Vazquez, 
arXiv:0805.4819v5 [hep-th]; Alejandro Vazquez, Hernando Quevedo and 
Alberto Sanchez, J.Gcom.Phys. 60: 1942-1949 (2010), arXiv:1101.3359vl 
[math-ph]; J.L. Alvarez, H. Quevedo and A. Sanchez, Phys.Rev.D 77: 
084004 (2008),arXiv:0801.2279vl [gr-qc]. 

[20] Hernando Quevedo, Alberto Sanchez, Sana Taj and Alejandro Vazquez, 
J.Phys.A A 45: 055211 (2012), arXiv:1101.4494vl [hep-th]; Safia 
Taj and Hernando Quevedo, arXiv:1104.3195vl [math-ph]; W. Janke, 
D. A. Johnston and R. Kcnna, J. Phys. A 43:425206 (2010), 
arXiv:1005.3392v2 [hep-th]; H. Quevedo, A. Sanchez, S. Taj and A. 
Vazquez, Gen.Rel.Grav. 43:1153-1165 (2011), arXiv:1010.5599vl [gr-qc]; 
Peng Chen, Int.J.Mod.Phys. A 26: 3091-3105 (2011), arXiv:1104.0546 [gr- 
qc]. 

[21] Deborah F. Jardim, Manuel E. Rodrigues and Stcphanc J. M. Houndjo, 
arXiv:1202.2830v2 [gr-qc]. 

[22] Jianyong Shen, Rong-Gcn Cai, Bin Wang and Ru-Keng Su, 
Int.J.Mod.Phys.A 22:11-27 (2007); arXiv:gr-qc/0512035. 

[23] Jan Aman, Ingemar Bengtsson and Narit Pidokrajt, Gen.Rel.Grav. 35:1733 
(2003); arXiv:gr-qc/0304015. 

[24] Hernando Quevedo, Gen.Rel.Grav. 40:971-984 (2008); arXiv:0704.3102v2 
[gr-qc]. 

[25] X. N. Wu, Phys. Rev. D 62: 124023 (2000). 

[26] Y. S. Myung, Y. W. Kim, and Y. J. Park, Phys. Lett. B 663, 342 (2008). 



22 



[27] Hernando Quevedo and Alberto Sanchez, JHEP 0809:034 (2008), 
arXiv:0805.3003v2 [hep-th]. 

[28] N. D. Birrell and P. C. W. Davies, Quantum fields in curved space, Cam- 
bridge University Press (1982). 

[29] G. W. Gibbons and S. Hawking, Phys. Rev. D 15: 2752-2756 (1977). 

[30] R. Hermann, Geometry, physics and systems (Marcel Dekker, New York, 
1973); G. Hernandez and E. A. Lacomba, Contact Riemannian geometry 
and thermodynamics, DifT. Gcom. and Appl. 8, 205 (1998). 

[31] G. Clement, J. C. Fabris and G. T. Marques, Phys. Lett. B 651: 54-57 
(2007). 

[32] Panagiota Kanti and John March-Russell, Phys.Rev.D 66: 024023 (2002); 
Wontac Kim and John J. Oh, J. Korean Phys.Soc. 52: 986 (2008); Kazuo 
Ghoroku and Arne L. Larscn, Phys.Lctt. B 328: 28-35 (1994). 

[33] S.P. Robinson and F. Wilczek, Phys. Rev. Lett. 95: 011303 (2005). 

[34] T. Jacobson and G. Kang, Class. Quant. Grav. lO:L201-L206 (1993); arXiv: 
gr-qc/9307002. 

[35] R. M. Wald, Genaral Relativity, Univ. of Chicago Press, 1984, Chicago. 
[36] L. H. Ford, arXiv: gr-qc/9707062. 

[37] Glauber Tadaiesky Marques and Manuel E. Rodrigucs, Eur. Phys. J. C 
72:1891 (2012), arXiv:1110.0079v2 [gr-qc]. 

[38] S.P. Robinson and F. Wilczek, Phys. Rev. Lett. 95,011303 (2005). 

[39] L. Smarr, Phys. Rev. Lett. 30: 71 (1973). 

[40] George Ruppeincr, Phys. Rev. D 75, 024037 (2007); Rev. Mod. Phys. 67: 
605-659 (1995), Erratum-ibid. 68: 313-313 (1996). 

[41] H. Quevedo, A. Sanchez, S. Taj and A. Vazquez, arXiv:1011.0122vl [gr-qc]. 

[42] R. Banerjee, S. Kumar Modak and D. Roychowdhury, arXiv:1106.3877v2 
[gr-qc]. 



23 



